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Abstract 

Using the technique of mean field theory applied to the lattice boundary 
Ising and tricritical Ising models we provide a qualitative description of their 
boundary phase diagrams. We will show this is in agreement with the known 
picture from boundary conformal field theory and we shall compare our work 
with that of Cappelli et al and show how their analysis is not in accordance 
with the physical picture. 



1 Introduction 



At criticality the lattice Ising and tricritical Ising models become conformally invariant 
and one can represent their behaviour using the M(3, 4) and M(4, 5) conformal field 
theories as is the case for all higher critical Ising models. On introduction of a boundary 
there have been found different boundary states with specific conditions that maintain 
conformal invariance. The simplest of these, and the ones we shall consider here, are 
the boundary states as found by John Cardy pQ known as the Cardy states. When 
considering the representative Ising and tricritical Ising lattice models, these Cardy states 
have physical interpretations in terms of the level occupancy for the boundary spins [2] . 
These being {±1} for the Ising model and {0, ±1} for the tricritical Ising model. These 
boundary states have been analyzed in the literature and it is known that there are 
relevant boundary fields present at some of the boundaries and on perturbation by one 
or more of these relevant fields a renormalization group flow will take one to another 
of the conformally invariant boundary conditions. For the Ising and tricritical Ising 
models under consideration here the details of the space of boundary RG flows has been 
found [3j and the flows express irreversibility as one would expect from the g-theorem 
[I], which states there exists a non increasing g- function on the space of boundary RG 
flows analogous to the c-theorem [5] in the bulk case. 

The boundary phase diagram for the Ising models under consideration has already 
been examined by Cappelli et al [6j in which they used a Landau Ginzburg type poly- 
nomial boundary potential. We however obtain a different picture based on a full mean 
field theory treatment and we will show when a Landau polynomial approximation is 
valid to the mean field equations and when it is not. The main present motivation for 
this is to develop techniques that can then be applied to the more complex case of the 
space of RG flows in the case of conformal defects. Furthermore as will be later eluci- 
dated not all of the flows in the boundary tricritical Ising model are integrable and our 
mean field description provides a qualitative picture of these flows. 

The layout of the paper is as follows, in the second section we will first outline the 
Ising space of boundary conformal flows and then we will derive the mean field phase 
diagram for both the bulk and boundary. In section three we will do exactly the same 
for the tricritical Ising model, in section four we will compare our picture with that of 
Cappelli et al and discuss when a Landau polynomial approximation is valid. In section 
five we will draw our conclusions. 

2 The Ising model 

2.1 Conformal boundary RG flows 

The Ising bulk conformal field theory comprises three symmetric holomorphic / antiholo- 
morphic representations of the Virasoro algebra and Cardy has shown that there are thus 
three conformal boundary states. Two of these states are stable in that they have no 
relevant boundary fields and these correspond to the boundary spins being completely 
fixed in either of the {±1} spin levels and these fixed boundary states are respectively 
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denoted (+) = (1, 1) and (— ) = (2, 1), where we have included the boundary states' Kac 
table index pairs. Finally one has the degenerate boundary state (d) = (1,2) where in 
the underlying lattice model fluctuations in the boundary spins take place equally across 
both {±1} levels. This degenerate boundary state has one relevant field 0i t 2 that the 
boundary magnetic field hf, couples to and switching on such a perturbation one induces 
an RG flow to either of the stable boundary states (±) with the fixed point values being 
hf, = ±oo respectively. Once at these fixed points as there are no relevant fields there are 
no more RG flows and the flows thus respect irreversibility as proposed by g-theorem. 
We now illustrate the simple space of flows where the dark dots indicate stable boundary 
states and the white dots represent unstable boundary states. 
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Figure 1: Conformal boundary RG flows for Ising model 
2.2 The mean field bulk model 

The square lattice Ising model we take as representative of the M(3,4) minimal model is 
composed of spins c?ij taking the values {±1}, i represents the horizontal lattice direction 
and j the vertical. We define the lattice action to be 

H = (-Ja i> j(a i+1 j + o-jj+i) - hdij) (2.1) 

Here J is taken to be a positive (ferromagnetic) nearest-neighbour coupling and h is an 
applied external magnetic field. We utilise the mean field approximation as described 
in [7J, this consists of first re-expressing the spin degrees of freedom a as deviations 
5a = (a — M) from their average value M resulting in a = 5a + M. Then the nearest 
neighbour interaction terms become a a' = aM' + a'M — MM' + 5a5a' and it is the 
last term, the product of fluctuations that one drops in the mean field approximation 
resulting in 

aa = aM + a'M - MM' (2.2) 



On entering this approximation into eq (2.1) one obtains the mean field action 



H mf = (-JvijMij + ^M id M id - hffijj (2.3) 

Here Mjj denotes the average magnetization at the lattice site and Mjj denotes 

the sum of nearest neighbour magnetizations of the lattice site ie. 

Mi,j = M ijj+1 + My-i + Mi+tj + Mi_ij (2.4) 
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From this we can deduce the mean field partition function 



Zmf = J>xp(-/3# m/ ) = Hexpi-^MijMij) 2 cosh(0( JM id + h)) (2.5) 

M i,j 

where the sum over {a} is taken over all lattice spin configurations and (3 = 1/T where 
T is temperature. The reason that the mean field partition function can be evaluated is 
that the sum over the spins decouple meaning that 

ex.p(-K n a n ) = J] ^ exp(-K n a n ) = 2 cosh(iT n ). (2.6) 

{cr n } n tr„=±l n 



From equation (2.5) we can obtain the mean field free energy 



fmf = -\ l og(Z mf ) = ^J2 M m M m " \ lQ g( 2 ™ S h(f3(JM id + h))). (2.7) 

Minimizing the free energy with respect to the magnetizations Mj results in the mean 
field consistency equations for the site magnetizations 

Mij = tanh(/3( JM id + h)) (2.8) 
2.2.1 Zero bulk fluctuation limits 

Before going on to study the nature of any phase transitions in the bulk we can use the 



previous magnetization equation (2.8) to deduce the regions of the phase diagram that 



correspond to all the bulk spins being completely aligned in either of the {±1} levels. 



Our first step is to take the previous bulk magnetization equation (2.8) and impose a 
flat magnetization profile as is expected in the bulk M(x, y) = fh. From this we can 
deduce the level occupancies p± that tell us the probability of finding a given spin in 
either of the {±} states from the relation p± = (1 ± in)/ 2 and this results in 

p± = ^(1 ± arctanh(/3(4Jm + h b ))) (2.9) 

On letting j3h tend to infinite limits (j3h — > ±oo) the fh term becomes insignificant and 
we can deduce the regions 

(Ph -»• ±oo) p± -»• 1. (2.10) 

These results are exact because far from any phase transitions mean field theory becomes 
exact. The reason for this is that going back to our original unapproximated Hamiltonian 



(2.1 ) and on taking the limits of h — > ±oo the unapproximated ha terms will completely 
dominate the nearest neighbour terms — J a a' that are finite on account of J being finite. 
One can deduce that if p± — >• 1 exactly in the limits (3h — >■ =too then any phase transition 
from the +1 level to the —1 level must take place for finite f3h. 
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2.2.2 Bulk phase transition 

Having established the stable regions where all the spins are aligned completely in either 
orientation we can now go on to study the nature of the phase transition between the 



levels. Going back to our original mean field equation (2.8) we want to consider the 
magnetization to be a smoothly varying field living on the plane and thus we now re 
express the site magnetizations in terms of a continuous function Mjj = M(ia, ja) where 
a is the lattice spacing. Next we Taylor expand the Mjj term to second order about 
the central coordinate (ia,ja). This then results in a differential continuum differential 



equation that approximates the bulk mean field equations ( |2.8| 

M = tanh(/3(4JM + a 2 JV 2 M + /i)). (2.11) 

This can be inverted to give the second order differential equation 

a 2 (3JV 2 M = arctanh(M) - 4/3 JM - f3h. (2.12) 

The rhs of the previous equation can be Taylor expanded for small M to give the second 
order equation as obtained in pQ 

(3a 2 JV 2 M = (l-4/3J)M + -M 3 - ph. (2.13) 

3 

In the bulk we impose a flat field profile with constant value denoted M(x, y) = fh Vx, y 
this in turn implies the second derivative term goes to zero and we can rewrite the 
previous equation as 

1 a 

= -0h+ (1 - 4/3J)m - -fh 3 . (2.14) 

3 

The condition for a critical point is a point of inflection on the curve of h against fh or 
dmh = d^h = 0; this condition is achieved at fh = 0, with h = and f3J = 1/4. The 



expansion of the rhs of equation (2.12) that resulted in equation (2.13) we will refer to in 
the rest of the paper as a Landau approximation. The reason being that these equations 
are the same as the ones which would have been obtained from the Landau-Ginzburg 
effective action 

S = f dxdy Q(VM(x, y)) 2 - phM(x, y) + \{l- A(3J)M(x, y) 2 + ^M(x, y) 4 

(2.15) 



The integrand in equation (2.15) would be obtained by Taylor expanding the free energy 
for small and slowly varying M(x, y), in such an expansion the coefficients of the different 
powers of M(x, y) have a functional dependence on the theory's parameters and it is the 
behaviour of these coefficients that dictates the nature of any critical points and phase 
structure around them. In the Landau theory, after making such an expansion one 
truncates the expansion at the smallest power of M whose coefficient stays positive 
under changes of the theory's parameters. The key point is that the order parameter M 
should be small for such an expansion to be valid. 
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2.2.3 Mean field boundary model 

Now we consider the lattice Ising model on the half plane and with the boundary running 
parallel to the vertical j direction and with the boundary's horizontal coordinate set at 
i = 0. The boundary conditions are set by a boundary magnetic field h b acting on the 
boundary spins ctoj and a ferromagnetic nearest neighbour boundary coupling J b that 
acts between the boundary spins, finally the boundary spins interact with the bulk spins 
via the bulk interaction J. These conditions are reflected in the following half plane 
Ising lattice action 

H = y~] (— Jaij(a i+1) j + (Tij+i) - hcfij) + (-Ja jaij - J b a 0> ja 0) j + i) - h b a j) . 

i>0,j j 

(2.16) 



We use the same mean field spin interaction approximation equation (2.2) as before for 
deriving the mean field partition function and further we assume that the only variation 
in magnetization will occur in the horizontal direction perpendicular to the boundary, 
meaning that Mij = = Mij+i V i,j and this results in a free energy per lattice 

row 

f mf =J h Ml + JM M 1 - 1 log(2 cosh(/3(2 J b M Q + JM 1 ) + h b )) (2.17) 
+ (jM l M l+1 - - log(2 cosh(/3( J(M,_! + 2M; + M i+l ) + h)))\ . 

Minimizing this with respect to the magnetization, df m f/dMi = Vz, we obtain the 
mean field magnetization equations. 

M = tanh(/?(2J fc M + JM{) + h b ) (2.18) 
Mj = tanh(/3(J(Mi„i + 2Mi + M l+ i) + /i)) (t > 0) (2.19) 

The following steps for obtaining a boundary equation for studying the boundary phase 
diagram are our own extension to the standard mean field techniques used in [7]. Firstly, 
as in the bulk case we again re-express the discrete magnetizations Mj in terms of a 
continuous function Mj = M(ia), where a is the lattice spacing and it is understood the 
magnetization function M(x) only depends on the horizontal plane coordinate x and 



thus the boundary magnetization value is simply M(0). The boundary equation (2.18) 
then becomes, on inversion and rearrangement, 

= arctanh(M(0)) - 2f3J b M(0) - fiJM(a) - (3h b . (2.20) 

It is this equation that we will use to obtain the boundary phase diagram, as it relates the 
boundary field h b with magnetization M(0) and any critical points are defined in terms 
of derivatives of h b with respect to M(0). The magnetization one lattice spacing from the 
boundary, M(a), in equation ( 2.20| ) will depend on the boundary magnetization M(0) 



via the bulk equations (2.19) and thus we need to obtain this dependence in order to end 
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up with an expression purely in terms of hf, and M(0). The strategy we shall employ 



will be the following: we first take the bulk equations (2.19) in their continuum form, 
Mi = M(ia) and we Taylor expand the magnetization Mi±\ = M(ia ± a) term about 
(ia) to second order. On inversion and rearrangement we obtain the same differential 
equation as we did in the bulk equation ( 2.12| ) 



a 2 f3Jd 2 M = arctanh(M) - 4/3 JM - (3h. (2.21) 

However, as there is by construction no variation in the vertical direction, we have in 
this case V 2 = d 2 where d will be used from now on to denote the derivative with respect 
to the horizontal direction. We take it that this differential equation exactly defines the 
form that M{x) will take and furthermore that it not only holds in the bulk (x > 0) as it 
came from the bulk magnetization equations but also at the boundary (x = 0) in order 
that we obtain a smoothly varying profile. As it is a second order equation we need to 
impose two boundary conditions in order to obtain the profile M(x). The first boundary 
condition will come from the fact that we impose a flat profile in the bulk limit x — > oo, 
the value of the magnetization in this bulk limit we denote lim x _ s>00 M(x) = Mm- This 
will entail that the second derivative term in our bulk profile differential equation will 
tend to zero in the bulk limit and we will be left with the bulk limit equation 

= arctanh(Moo) - I U.U X - (5h (2.22) 



This is the same as the bulk theory equation (2.14) obtained in the previous section and 
thus we find the same critical parameters of h = and f3J = 1/4 with the result Moo = 0. 
As we are looking to hold the bulk critical while varying the boundary parameters, we 



shall insert these critical bulk values into our profile differential equation (2.21) in order 
that we obtain a critical bulk profile equation 

a 2 d 2 M = 4(arctanh(M) — M). (2.23) 

As previously mentioned, one of the boundary conditions will be the critical bulk limit 
Mqo = and the second will be the value of the boundary magnetization M(0). We can 
integrate this second order equation once including the bulk limit boundary condition 
Mqo = to obtain a first order critical profile equation with only one boundary condition 
the boundary magnetization M(0) to be set. We use the result 

1 r M 

-(dM) 2 = / d 2 M dM (2.24) 

2 Jo 

to find 

adM(x) = ±2 v / 2Marctanh(M) + log(l - M 2 ) - M 2 . (2.25) 
As shown in Cappelli et al [6] the correct choice of sign in equation (2.25) is positive 



if the value of the magnetization is negative as the magnetization should monotonically 
increase to the bulk limit value of zero and vice versa. So this first order equation defines 
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the form of M(x) given a boundary value of M(0) and this means that the M(a) term 
in the boundary equation can be expressed in terms of M(0) in the following way. Given 
that we have the values of the first and second derivatives at a point we can use them 
together to obtain all higher order derivatives using the chain rule 

8 n M = d M (d n ~ 1 M)dM. (2.26) 

And with knowledge of all the derivatives we can write down a series expansion for 
M(a) about the boundary resulting in an expression purely in terms of the boundary 
magnetization M(0) in the form of an expansion B e {M (0)) 

M(a) = B e (M(0)) = £ a "^ (0) . (2.27) 

n=0 

Even though we can never consider the whole expansion we will find that at critical 
points where M(0) = we will only need to consider the behaviour of the first two 
terms (the behaviour of the rest follows by induction) and when M(0) 7^ as occurs 
when considering the tricritical model we can suitably truncate the expansion to a degree 
that does not appreciably affect any numerical answers. With this expansion and on 
letting M(0) = m our original boundary equation becomes 

= arctanh(m) - 2(3J b m - ^B e (m) - /3h b . (2.28) 

This is an expression purely between h b and m that we can use to study the boundary 
phase diagram. Before we go on to obtain the boundary phase diagram we will point out 
that we cannot use the above equation to plot m vs h b across the whole range — 1 < m < 1 
as the expansion B e (m) only converges for the range \m\ < 0.6. To overcome this we 
will obtain a numerical form of the function M(a) = B e {m) denoted M(a) = B n {m) 
and will use it when making plots of the full mean field equations across the whole range 
of boundary magnetization. The way we obtain this function is to numerically solve 



our critical first order profile equation (2.25) a large number of times each time with a 



different value of boundary condition m taken evenly from the full range — 1 < m < 1. 
Then for each of these different m dependent numerical solutions we read off the corre- 
sponding value of M(a) and tabulate the results in a list of M(a) vs m and then we use 
Mathematica to define a interpolative numerical function M(a) = B n (m) from said table. 
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Figure 2: Plot of numerical function M{a) = B n (m) that gives the magnetization a 
lattice spacing from the boundary in terms of the magnetization at the boundary 



With this we obtain a numerical boundary equation that can be used to plot the full 
mean field behaviour of m(h b ) 

= arctanh(m) - 2/3J b m - ^B n (m) - (3h b . (2.29) 



However it is the expansion form of the equation (2.28) that we will use to deduce the 
existence and nature of any critical points as will follow. First though we examine the 
regions of the phase diagram where there are no boundary fluctuations. 



2.2.4 Zero boundary fluctuation limits 

In analogy with the bulk there are limits of the boundary phase diagrams that correspond 
to fixing all the boundary spins either of the {±1} levels. In this case we need the 
boundary occupancy probabilities p± which we obtain from the same procedure as in 
the bulk to be 

P± = -(1 ± arctanh(/3(ET + h b ))) (2.30) 

where Kq = 2J b M(0) + JM{a) and this term is always small and finite and thus as in 
the bulk it is completely dominated by any large value of h b and we find the regions 

(f3h b -»• ±oo) p± -»• 1 (2.31) 

and we can identify these respective limits j3h b — > ±oo with the stable Cardy states (±) 
that correspond to aligning all the boundary spins in one orientation. 



2.2.5 Boundary phase transition 

As in the bulk, we can expect a phase transition at h b = m = that separates the stable 
limits considered before. As a starting point we take the boundary equation (2.28) 
obtained before 

= arctanh(m) — 2(3 J b m — -B e (m) — f3h b . 
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We cannot expand the rhs of the equation (2.28) in powers of m as is normally the case 



due to the behaviour of the odd terms in B e {m) coming from taking different signs of 
the square root in our first derivative term. We can, though expand, in powers involving 
\m\ resulting in 



= -Bh h + (- - 2BJ b )m -\ -=m \m\ + 0(m 3 \m\). 

4 2^/6 



(2.32) 



We obtain a phase transition for hy, = m = when the linear term in equation (2.32) 
vanishes at ftj^ = 3/8 and this entails d m hb = as is required for a critical point. 
However, normally, one requires d^hb = as well, but the second derivative is not 
defined at m = in the same way that the first derivative of \m\ is not defined at m = 0, 
though graphically it at least looks like we have a point of inflection. This critical point 
we identify with the (d) boundary state from the known conformal picture. In figure 
([3]) we plot the phase transition using both the numerical form of the full mean field 
equation and also the polynomial Landau approximation where it should be noted that 
we have to set the bulk interaction J and we take it to be J = 1 in the all the plots in 
this paper. 
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Figure 3: Plots of m(hb) using (a) the full mean field equations and (b) the Landau 
approximated version of the mean field equations for a smaller range of m as the Landau 
approximation only holds for small values of m 

The full mean field plot qualitatively illustrates the picture from the conformal field 
theory with RG flows taking one from the (d) state at hf, = to either of the stable states 
(±) with m = ±1. From now on in we shall only plot phase transitions using the full 
mean field equations as the Landau approximation plots just elucidate the behaviour 
of the full mean field plots for small values of magnetization ie they don't add any 
information. 



2.2.6 Superpotential 

It is sometimes useful to consider the boundary equation the master boundary equation 
(2.28) as coming from the minimization of a function of the boundary magnetization 
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m. This was first considered by Cappelli et al in [6] where they called this function 
the '"superpotentiaF". Accordingly we define our superpotential W(m) such that the 



condition d m W = gives rise to our master boundary equation (2.28) . This is only 
possible exactly when taking the Landau polynomial approximation of the equations as 
the terms in B e {m) are not amenable to exact integration. However we can numerically 
obtain a superpotential W(m) of the full mean field equation as B n (m) can be integrated 
numerically, this allows us to make plots of the superpotential across the whole range 
— 1 < m < 1 and also when a Landau type expansion is not possible as occurs at 
some critical points of the tricritical model. Taking the Landau form of our boundary 
equations and integrating we obtain 

W Lan (m) = -ph h m + (| - (3J b )m 2 + -^=m 2 \m\ (2.33) 

if we tune to its critical value of 3/8 then this expression is equivalent to the super- 
potential expression obtained by Cappelli 

W cap (m) = aom + -m 2 \m\ (2.34) 

with — (3hb = ciq. As stated before we can also obtain a mean field superpotential 
W m f(m) taken by integrating the full mean field boundary equation with the numerical 
B n {m) integrated numerically, this results in 

1 If 

W m f(m) = marctanh(m) + - log(l — m 2 ) — fiJ^vr? — — / B n (m)dm — f5h[,m (2.35) 

In figure Q we plot the superpotential in both full mean field form W m f(m) and in the 
polynomial approximation. 



W(m) W(m) 




(a) -i o i (b) -0.05 o 0.05 

Figure 4: Plots of superpotential W(m) for large (a) and small (b) m, the dashed plot 
corresponds to the Landau polynomial approximation and the full plot corresponds to 
the full mean field superpotential, note that in figure (b) plots are coincident near origin. 



as was the case for the phase transition plots we shall only plot the full mean field 
superpotential from now on. It should be noted that in subfigure (a) figure @ the 
vertical axis increment of 0.238 is chosen as W m f(H) = 0.238 3sf. Outside of the range 
— 1 < m < 1 one finds that W m f(m) becomes complex and unphysical and therefore the 
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mean field picture naturally retains the magnetization in this range, however the Landau 



approximated superpotential (2.33) can be plotted across all values of m and does not 
have this retantion property. In fact all the mean field superpotentials considered in this 
paper have this retention property. An advantage of this is that it allows one to deduce 



the level probabilities (2.30) for large m whereas in the Landau picture where there is 
no such retention one loses this information. In figure ([5]) we illustrate the boundary 
Ising phase diagram using our full mean field superpotential. 




Figure 5: Mean field Ising boundary phase diagram 



3 The tricritical Ising model 

3.1 Conformal RG boundary flows 

The tricritical bulk theory is described by the M(4, 5) minimal model and has six holo- 
morphic/antiholomorphic representations of the Virasoro algebra and thus there are six 
corresponding Cardy states when one introduces a boundary. As in the case of the Ising 
model these can be related to the conditions on the boundary spins of the represen- 
tative tricritical ising lattice model where the spins can take the values {±1,0} which 
correspond to spin up/down and spin zero or no spin respectively. There are stable 
boundary states that have no relevant boundary fields, these correspond to fixing all 
the boundary spins in a given orientation. They are denoted as follows where we also 
include their corresponding Kac table indexes: (+) = (1,1), (— ) = (3,1), (0) = (2,1). 
We also have two mixed boundary states (0+) = (1,2) and (0—) = (1,3) that corre- 
spond to allowing the boundary spins to only fluctuate between the levels {0, +1} and 
{0, —1} respectively. These states have one relevant boundary field and perturbing 
with this field will induce an RG flow to one of the given stable states as will be shown 
in the next figure. The last of the Cardy boundary states is the degenerate boundary 
state (d) that corresponds in the lattice model to letting the boundary spins fluctuate 
across all levels. It has two relevant boundary fields 0x,2 an d 4>i,3, on perturbations by 
these fields one obtains an RG flow to the other Cardy states as is elucidated in figure (T6J) 
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Figure 6: Space of boundary RG flows from Ref. [3], finely dashed lines represent 1st 
order transitions, coarsely dashed represent non integrable flows 

As one can see from figure Q the mixed flows taking one to either of the (±0) states 
are not integrable and hence our following mean field boundary phase diagram will give a 
qualititive picture of them. One can also consider a superposition of the stable boundary 
states (+) © (— ) as first show in [8]. This superposition has one relevant boundary field 
given by the second identity field, on perturbation by this field one obtains a first order 
transition (a discontinuous jump in magnetization) to either of the stable states (+) 
or (— ). In the next subsection we shall apply mean field theory to the tricritical Ising 
lattice model in order to obtain a qualitative picture of the boundary phase diagram and 
then we shall see how this compares with the conformal picture just described. 



3.2 The mean field bulk model 

Now we consider the tricritical Ising model where we allow the spins a to take the value 
{0} as well as {±1}, alternatively the {0} state can be viewed as a spin vacancy. The 
square lattice action we shall take is 

H = ^2 (-J(Tij(<Ti+lj + (Tij+l) - haij + Aofj) . (3.1) 

y 

The parameter A is known as the vacancy density and increasing it in the positive 
direction increases the number of vacancies and increasing it in the negative direction 
decreases the number of vacancies. The parameters J and h are as before a ferromagnetic 
nearest neighbour coupling and an applied magnetic field. On applying the same mean 



field substitution (2.2) of the nearest neighbour terms in the action we obtain the mean 
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field action 



H mf = Yl (- J °i,i M i,i + \ M ij M H - h °i,3 + . (3.2) 

From this we obtain the mean field partition function 

Z m f = X\eM~M^M hj ) (1 + 2y cosh(/3( JM itj + h))) (3.3) 

where y = e~^^ and from this follows the free energy 

fmf = ^Yl M ^Mi,j - - log(l + 2y cosh(/3(JM M - + h))). (3.4) 

Minimizing the free energy with respect to the magnetizations Mjj results in the mean 
field consistency equations for the magnetization 

Mi = 2 y sin HP(JM ijL + h)) 
hJ 1 + 2y cosh(/3( JiWj j + h)) 

before using this equation to study the bulk phase transitions we will use it to study 
regions of the phase diagram where all the spins completely occupy one of the {0, ±1} 
levels ie where there are zero fluctuations. 

3.2.1 Zero bulk fluctuation regions 

In order to deduce the regions we need the level occupancies p± and po obtainable from 
the relations po = 1 — (c 2 ) and p± = ((c 2 ) ± (<r))/2. We can obtain the expectation 
value of (<r 2 ) by differentiating f m f wrt to A and we obtain it as the same expression for 
the magnetization expect with the sinh turning to a cosh. In the bulk theory we impose 
a flat magnetization profile and this means that Mij = fh Vi, j where fh is the average 
bulk magnetization. On inserting this into our magnetization and (c 2 ) equations and 
following the previously mentioned relations we obtain the level occupancies 

e [3(4Jm-(A-h)) 

(3.6) 



1 + e /3(4Jm-(A- 


-h)) _|_ e /3(-4Jm- 


-(A+h)) 


e P(- 


4Jm-(A+ft)) 




l + e /3(4Jm-(A- 


-h)) _|_ e /3(-4Jm- 


-(A+h)) 




1 





(3.7) 



P0 1 + e f3(4Jm-(A-h)) + e l3(-AJrh-(A+h)) ' ( 3 " 8 ) 
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Using these expressions we deduce the zero fluctuation regions as 



/3(A — /i) — y oo and /3(A + h) -> oo =>- p -> 1 (3.9) 
/3(A - /i) -)■ -co and /3/i -)■ oo =>- p + ->■ 1 (3.10) 
/3(A + /i) -> -co and -> -co => p_ -> 1 (3.11) 

We will analyze these regions in more detail in the boundary case where the analysis is 
essentially the same as in the bulk case. 



3.2.2 Bulk phase transitions 

Now we will examine the phase transitions present in the bulk model and for this we 



need to consider our previous magnetization equation (3.5). This can be inverted, the 
Mij Taylor expanded and the Mi j can be recast as M(ia,ja) as before resulting in 

a 2 (3 JV 2 M = F(M, y) - 4/3 JM - (3h (3.12) 

where the function F(M, y) is defined as follows 



+ (-3) 



The magnetization profile M = M(x,y) must obey the differential equation (3.12) in 
the bulk, however as there is no variation in the bulk we take the profile to be again 
constant M(x, y) = rh Vx, y this then gives 

0=F(fh,y)-4pjfh- f3h. (3.14) 

This can be expanded in a Landau approximation for small m giving 



3 



+ (3- 18i ,,3V)(l + 2,)3 ^ + 0(rn6) 

This is critical if dmh = d^h = 0, which is equivalent to the linear and constant terms 
going to zero in the above ie h = and (3 J = (1 + 2y)/8y. Furthermore it is tricritical 
if <5%h = which corresponds to the cubic term vanishing as well which happens for 

h = 0, y= i PJ= | (3.16) 

Thus we have a tricritical point at these coordinates and from here increasing y, which 
corresponds to decreasing A takes one along a line of critical points and in the limit 
y — > oo A — > -co we recover the Ising bulk mean field equation and in turn the critical 
Ising phase transition. Furthermore we can find two more lines of critical points adjoining 
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the tricritical point, in fact it is the joining of the three lines of critical points that gives 
the tricritical point its name. These two lines of critical points occur roughly along the 
lines (h = ±A) and we shall only consider the limits (h — > ±00, A — > 00 =>■ y — > 0) as it 
will be these limits that are of interest when we will consider the boundary theory. On 
taking these limits our bulk equation takes the form 

= 2arctanh(2m T 1) — 4/3Jm - (3(h T A). (3.17) 

These equations cannot be expanded for small fa as the arctanh term diverges at m = 0. 
However we can still solve it exactly and one finds critical points when fa = ±1/2, 
(3 J = 1 and (3{h =F A) = =f2. Having found the critical values of fa = ±1/2 we can 
Landau expand around these critical points by defining m± = fa =F 1/2 

16 , 

= -p(h T A T 2J) + (4 - 4/3J)m± + —m\ (3.18) 
3.3 The mean field boundary model 

As in the case of the Ising model we consider the tricritical model defined on the infinite 
half plane with lattice action 

H =Y1 {- Ja i,j( a i+hj + a i,j+l) ~ h(J iJ + Aa i,j) ( 3 - 19 ) 
i>0,j 

+ ^2 (~ Jcr o,jO-i,j - Jbcr ja 0tj+1 - h b a ,j + A^oj) 
3 

here the boundary parameters are denoted with subscript b. We assume the magnetiza- 
tion only varies perpendicular to the boundary and on applying the mean field approxi- 
mation we obtain the free energy 

f mf = J h Ml + JMqMi - 1 log(l + 2y b cosh (/?( 2 J b M + JMi) + h b )) (3.20) 

+ (jM l M l+1 - - log(l + 2ycosh(/3( J(Mi_! + 2M t + M i+1 ) + /i)))") . 

On minimizing the free energy with respect to the magnetization we obtain mean field 
consistency equations for the boundary and bulk magnetization. 

_ 2y b smh(P(2J b M + JM 1 + h b )) 

l + 2y b cosh((3 (2 J b M + JM 1 + h b )) 1 ' j 

= 2y sinh(/3( J(M,_x + 2M, + M m ) + fe)) 

1 1 + 2y cosh(/5( J(Mi_i + 2M, + M m ) + fc)) 1 ^ j 1 ' 

For the bulk equations we impose the tricritical bulk parameters h = 0, y = 1/4 and 
(3J = 3/4 as we are looking to model the tricritical Ising conformal field theory in the 
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bulk. Then we invert them and on going to a continuum description Mj = M(ia) we 
obtain 



a 2 d 2 M = - log 



2M 
1 - M 



+ 



1 + 3M 2 
M 2 — 1 



4M. 



(3.23) 



This is equivalent to the bulk differential equation (3.12) obtained in the last section 
tuned to the tricritical point and with V 2 = d 2 as we assume no profile variation in the 
vertical direction. The boundary equation on inversion and in the continuum picture 
becomes 



= F(M(0), y b ) - 2/3 J 6 M(0) - pjM{a) - (3h b . (3.24) 

As with the Ising case we need to convert the above boundary equation into an expression 
containing only the boundary magnetization that we relabel M(0) = m. Thus we need to 
express the M{a) term as a function of m, M(a) = B{m). We shall define this function 
in two different ways in exactly the same manner as was with the Ising model, the first 
as a series expansion B e (m) and secondly as a numerical function B n {m). We shall use 
the series expansion to analyze the nature of any critical points but as the expansion 
does not converge for \m\ < 0.55 we shall use the numerical function in order to produce 
any plots of m(h b ) or mean field superpotential W m f(m) in order that we can plot them 
across the whole range of m. Inserting the expansion form B e {m) back into the equation 
(331) 



F(m, y b ) - 2/3 J b m - - ■ B e {m) ~Ph b = (3.25) 

we obtain a master boundary equation relating the boundary magnetization M(0) = m 
and field h b that we can use to deduce the existence and nature of any critical points. 
Though first we shall investigate the parameter limits we expect to correspond to the sta- 
ble conformal boundary conditions ie limits where there is no fluctuation in the boundary 
spins' occupancy. 



3.3.1 Zero boundary fluctuation limits 

In order to investigate the behaviour at the boundary parameter limits that one would 
expect to correspond to the stable conformal boundary states (0) and (±) we need the 
mean field occupancies po,p+ and p_ that tell us the probability that the boundary spin 
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will occupy the or ±1 levels. We report them as 



e /3{K -(A b -h b )) 



(3.26) 



P+ = 



1 + e /3(K Q -(A b -h b )) + e /3(-K -(A b +h b )) 



e P(-K -(A b +h b )) 



(3.27) 



P- = 



1 + e /3(K -(A b -h b )) + e 0{-K o -(A b +h b )) 



1 



(3.28) 



Po = 



1 + e P(K —(A b —h b )) + e P(-K -(A b +h b )) 



where Kq = 2Jt,Mo + JM\ and this term is always finite and thus will be dominated 
by infinite limits of ± hb The conformal boundary RG flows end in infinite limits 
of the boundary parameters hb and A& and thus we look for the mean field behaviour 
at the end of curves hb(Ab) in the hb vs plane. There are two cases to consider, 
the first is when A5 < and when A& — > —00 we get po — > for any curve hb(Ab). 
Furthermore if limA () -»-oo/i&(A&) = ±00 we get p± — > 1 regardless of the how the curve 
hb(Ab) achieves this limit. In order that we have a phase transition between the ±1 
boundary states we need neither p± — > 1; the only curves that make this possible are 
ones such that limA^-oo^A^) = constant, ie curves that are asymptotically parallel 
to the line hb = 0. In this respect the line /if, = for A& < can be seen as a boundary 
between the p± = 1 stable limits. 

Now we consider the case such that A5 > 0. When A5 — > 00, studying the occupation 
probabilities we see that for curves hb(Ab) st limA^oo ( A& — hb(Ab)) = —00 ie for curves 
whose asymptotic gradients are greater than 1 we end up with p + — > 1. Similarly for 
curves st limA 6 ->oo(Ab + hb(Ab)) = — 00 ie ones whose asymptotic gradient is less than 
1 we get p_ — > 1 .For curves liniA^oo ( Aj, =F hb(Ab)) = 00 ie ones whose asymptotic 
gradients are between and not including plus or minus one we get po — > 1. From this we 
can draw the conclusion that any phase transition between the level and the 1 level must 
occur at the endpoint of a curve that is parallel to the line hb = A5. Similarly for a phase 
transition between the and —1 levels must occur at the end of a curve asymptotically 
parallel to the line hb = —Aft. Our next figure illustrates how the boundary parameter 
plane is divided into the respective regions. 
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(+) 

h b =_0_(A b ^<0) 
(-) 



h b =A b (A b >0) 



(0) 



-► 



h b =-A b (A b >0) 



Figure 7: Regions of no boundary fluctuations 



These results can be taken to be exact because mean field theory becomes exact in 
the regions of the phase diagram far from any critical points. As in the Ising case if 
we go back to our original lattice action (3.19) and take these stable limits, the terms 
—hb&ij and A&cr? • dominate completely the nearest neighbour terms. 



3.3.2 Boundary phase transitions 

In correspondence with the (0±) Cardy boundary states we anticipate critical points in 
the limits of (hb — > ±oo, Aj, — > oo) along a curve parallel to the line hb = ±A?,. We 
shall denote these lines hb = ±A& + u>± (note that here we revert to using instead of 
Ub = e~^ Ai >). Taking these respective limits (hb — > ±oo, — > oo) our master boundary 
equation takes the form 

3 

2arctanh(2m =F 1) - 2/3 J b m - -B e (m) - f3u± = 0. (3.29) 

As in the bulk we are unable to expand this equation around m = as the arctanh 
term diverges there, in the bulk we obtained an exact value for the critical value of 
fh = 1/2 and we could then retrospectivley expand around that. In the boundary case 
as we shall see at best we are able only able to obtain an approximate numerical value 
for the critical value of the boundary magnetization and as a result we cannot expand 
around it. Thus we use the point of inflection condition of d m hb = d^hb = on the 
full mean field equation as the condition for criticality. Obviously we cannot consider 
all the terms in the expansion B e (m) and so we truncate it to the first twenty terms 
as including more makes no difference to the numerical answers obtained. In doing so 
we obtain critical points at these limits with the parameters taking the critical values 
m = ±0.477, f3Jb = 1.76, w-t = =f2.77 (3sf). In figure ^ we plot the behaviour of m(ui±) 
and W(m) at the respective (0±) limits just explored. 
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(a) -tot (b) 

Figure 8: Plots of (a)m(w+) at (0+) limit (b) m(u;_) at (0—) limit 



These plots indicate the fact that u± are the parameters defining the boundary 
conditions once we have taken the limits A& — > oo, hb/Af, — > ±1 that correspond to the 
respective (0±) conformal boundary states. 

The conformal space of flows would suggest that a first order transition should be 
present in the limit (A/, — > — oo, h b = 0). From our investigation into the stable limits 
we concluded that there would be a phase transition from the {+} level to the {— } level 
in this limit of (A^ — > — oo, h b = 0). In this limit our master boundary equation (3.25) 
becomes 



arctanh(m) — 2(3 J^m B e (m) — fih b = 



(3.30) 



this can be expanded for small m resulting in the Landau expansion 



(3.31) 



Here one obtains a critical phase transition when h b = and /3Jf, = 1/8. However as 
previously stated results from conformal field theory tell us that there should be a first 
order transition (discontinuous jump in magnetization) present here. In order for a first 
order transition to occur in our mean field theory treatment we would require f3 J b > 1/8 
and one gets three solutions to the mean field equation, one being m = and the others 
m = ±a where < a < 1. Normally one would not consider the zero solution as physical 
as it sits atop a maximum of the potential, however spontaneous magnetization is not 
possible at a one dimensional boundary and thus in our modeling we will take m = 
to be the physical solution for /if, = and when a small hf, is switched on the boundary 
magnetization jumps to either m = ±q depending on the sign of the perturbation. In 
this sense we effectively model the exact picture coming from conformal field theory. It 
should be noted however that in our mean field picture if one wanted a maximal first 
order transition ie a finite jump to either m = ±1 then fij b would have to tend to 
infinity or conversely the temperature would have to be zero. In figure ^ we produce 
full plots of both m(h b ) and W(m) setting f3J b = 2/5, quite an arbitrary value but one 
that illustrates the possibility of a first order transition. 
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(b) 

Figure 9: Plots of m{hb) when Aj, — > oo (a) for /3J& = 2/5, the black dot indicates the 
zero solution we take as being physical in our interpretation (b)for j3J{, = 1/8 

We have considered the phase transitions occurring at the limits of lines parallel to 
hb = iAft, and at most one could obtain a second order critical point along said lines. 
We expect a phase transition at roughly the intersection of these three lines in analogy 
with the bulk tricritical point. Expanding the lhs of our master boundary equation 
(3.25) for small m we obtain 

Ph + 







1 



1, 



+ 



2y b W + ^m+l 
(3-18y b + 32y fc 2 )(l + 2y fe ) 3 



(4y b -!)(! + 2y b ) ; 
48y b 3 



+ 



3 3 



m 



(3.32) 



— (56 + 5VT5) ) m 5 + 0(m 7 ) 
560 / 



1280 % 5 

this is critical if f3 J b = (2 + y&) /8yb and hb = and tricritical if 

Vb = Ubtr 



1528200 + 388800V15 - 16200V 17361 + 4528V15 



240 + 108V15 



+ 6 ■ 53 283 + 72V15 + 3\/ 17631 + 4528V15 



60 



=0.216 (3sf) (3.33) 

This tricritcal phase transition we identify with the conformal (d) state. We now illus- 
trate the tricritcal phase transition in figure ( 10 ) 
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Figure 10: (a)m(hb) through tricritical point (hb = 0, A& = — 4/3 log(yfc cr ^) = 2.04(3sf)) 
indicated by a black spot (b) (c 2 ) (A;,) through tricritcal point again indicated by a 
black spot 



We have included a plot of (cr 2 ) versus A& for /i& = 0. This we hope illustrates the 
qualititive nature of the conformal boundary flows from the (d) state to either of the 
(+) © (— ) and (0) states. With (a 2 ) = 1 corresponding to the state (+) © (— ) and 
(V 2 ) = corresponding to the (0) state. In figure (11) we summarize the boundary 
phase diagram using our full mean field superpotential W m f(m) obtained by integrating 
the numerical form of our boundary equation as in the Ising case. 
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Figure 11: Boundary tricritical Ising model phase diagram 



In the center we have the tricritical point we identify with the conformal boundary 
(d) state, which similarly to the bulk occurs at the intersection of three lines of critical 
points. We only explicitly worked out the case for the infinite limits of said critical lines 
but on investigation one finds that there are critical points along the whole of the lines. 
The limit of the two of these critical lines with A;, — > oo we identify with the (±0) 
conformal states. The — > — oo critical line we plot with > 1/8 in order to give a 
first order transition as is expected from conformal field theory and we identify this limit 
with the (+) © (— ) conformal boundary state. Any other curves ending in infinite limits 
of the boundary parameters we have shown will end up at stable boundary conditions 
that we identify with the (0), (±) conformal boundary states. The smaller plots in the 
figure relate to perturbing the boundary parameters slightly about the phase transition 
in question. 
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4 Comparison with work of Cappelli et al 



In the case of the Ising model considered in section 2 we saw that Cappelli's superpo- 
tential agreed with our Landau approximation to the boundary mean field equations 
if we tuned fiJ b to it's critical value. However such an approximation is only valid for 
small values of the order parameter m, in the regions h b — > ±00 we get m — > ±1 of 
the phase diagram where such a Landau approximation of the full mean field equations 
would not hold and furthermore in these regions mean field equations become exact. If 
one does take the Landau approximation to hold in all regions one would get a picture 
where m — > ±00 and it would be these limits that Cappelli identifies with the confor- 
mal boundary states (±) arguing that in the conformal picture any nonzero boundary 
magnetization would get renormalized to infinity and thus it is natural to identify these 
limits. 

In the case of the boundary tricritical Ising model Cappelli put forward the following 
polynomial potential W ca p(m) 

2 4 4 

rrr I \ 1 771 171 , f 2 171 /A -t \ 

W cap (m)=bm-a— + — = -h b m + d b rn +— (4.1) 

where in the second equality we relabel his parameters b = — h b and 5 b = —a/2 for 
direct comparison with our work. He identifies certain points and limits in the boundary 
parameter space h b vsS b with the the Cardy boundary states as we have done. In his 
analysis he states that the boundary value m should be given as the stationary points 
of his superpotential, the stationary equation d m W cap = is the following cubic 

= -h b + 25 b m + m 3 (4.2) 

one can differentiate again to find the equation that defines a point of inflection for 
W cap , d^Wcap = 0, this results in m = ±(— 25 b /3) 1 ^ 2 at these points of inflection and 
on substitution of this value of m back into the stationary equation we obtain the 
relation between h b and 8 b as is necessary for an inflection point. One obtains h b = 
±2(-2(5 b /3) 3 / 2 , for S b < these wings define boundaries in the parameter space, to the 
right of the wings there is only one real solution and two complex conjugate solutions to 
the cubic stationary equation and the real solution will be a minimum. On the wings and 
to the left of them there will be three real solutions with different degrees of degeneracy. 
To the left of the wings there are three distinct real solutions, on the line h b = 0, 6 b < 
one of the solutions will be m = and the other two will be symmetric about this point 
m = ±a. Please refer to the next figure for the following discussion, (the wings are 
dotted lines). 
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Figure 12: Tricritical Ising model phase diagram of Cappelli et al 



On the wings themselves (not including origin hb = 5b = 0) two of the distinct 
solutions present just to the left of the wing will have coincided on going onto the wing 
and the coincident solution represents the point of inflection. The third real solution on 
the wing will be a minimum. At the origin hb = = there will be three coincident 
solutions of m = 0. Cappelli uses the degeneracy of the stationary points as the guiding 
light in his analysis as to what points or limits to identify with which Cardy states. The 
reason being is that there can be a nice match between the number of relevant fields 
at the given Cardy states and (n — 1) where n is the degeneracy of the points/limits 
Cappelli chooses. Accordingly as at the origin n — 1 = 2 and as there are 2 relevant fields 
at the (d) boundary state he makes the identification of the origin hb = 5b = with the 
free boundary state. As n — 1 = 1 along the wings and there is one relevant field at the 
(±0) states he identifies them with the limits (hb — > ±oo, 5b —> — oo) along the wings. 
To the right of the wings n — 1 = and he thus identifies certain infinite limits of hb and 
5b here with the the stable boundary states (±), (0). He does not however specify along 
which curves one should take these limits. Only that any limit (hb — > ±00,^ — > 00) 
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should be identified with the respective (±) states and that the limit (hb = 0, 5b — > oo) 
should be identified with the (0) state. For the Cardy state superposition (+) © (— ) 
Cappelli identifies the limit (hb = 0, 5b — > — oo) however the match between n — 1 = 2 
here and the number of relevant boundary fields does not hold, possibly one could say 
identify solutions of the form m = ±a. 

We think there are a number of issues with this picture of the boundary phase dia- 
gram and we will now compare Cappelli's picture with our full mean field picture and 
argue that ours gives the correct physical picture. Starting with the phase transition cor- 
responding to the (d) state in our picture we obtain the polynomial expansion WLanim) 
by integrating the our Landau expanded boundary equation ( |3.32[ ) 

WLm _ _ Phbm + i ( _L _ 2M + i )ra . + y*-W + *r + (4.3) 



1 /(3-18y 6 + 32y b 2 )(l + 2y^ 27 



+ - — -(56 + 5V15) m° 

6 V 1280y£ 224CT 'J 

the coefficients in this expansion have a complicated functional dependence on A&, for 
a more direct comparison we can expand the coefficients around the tricritical values of 
the parameters 

(3J b =yfl±l+ e p Ab = - log{ybtri) + S (4.4) 

to result in 

W Lan = -h b m + (<5 6 - e)m 2 - 0.474<W + (0.293 + 1.045 b )m 6 (4.5) 

the decimals in the coefficients can be given as very long exact answers but we give 
present them as decimals given to (3sf). Now our term 5b is directly comparable to 
Cappelli's and we see that if we tune e = then the coefficients of our m 2 terms become 
the same, but again there is further coupling at the higher order terms. This entails 
that we obtain a tricritical point in analogy to the bulk picture when we tune 5b = 
whereas Cappelli's analysis predicts a critical point. It is important to state that by 
construction our potential WLan is only valid for values of the boundary parameters 
close to the tricritical point. 

The next issue is with the (±0) limits, in Cappelli's picture he has the unphysical 
behaviour in that he has the (0+) state occurring for a negative hb, whereas the (+) 
state occurs for positive hb, meaning that the magnetization field changes the nature 
of its alignment with the external field at different regions in the phase diagram. In 
our picture we obtain consistent alignment between the magnetization and external field 
throughout the phase diagram. Furthermore we found that we could n't even obtain a 
polynomial Landau expansion and had to rely on the full mean field equations at these 
phase transitions, even in the bulk case where we could obtain a polynomial expansion 
we had to take the limits (hb —> ±oo, — > oo) first in the full mean field equations and 
then make the expansion. There was no way we could obtain a polynomial expansion 
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that held in all regions of the phase diagram. This is to be expected as the crux of 
the Landau approximation is to expand around a small value of the order parameter 
and in the phase diagram of the bulk or boundary models there are regions where the 
order parameter is not small. The final issue with Cappelli's picture is that our analysis 
showed exactly that any phase transition (±0) that occurred at the infinite endpoints 
of curves in the h b vsA b plane had to occur at the endpoints of lines asymptotically 
parallel to the lines h b = ±A&. Cappelli says the (±0) occur at the endpoints of the 
wines h b = ±2(-2(5 fe /3) 3 / 2 , unless there is some very particular coupling between the 
vacancy density as defined in the lattice Hamiltonian and the quadratic coefficient 5 b 
in Cappelli's picture it is hard to see how the wings would n't end in stable limits. 
Unfortunately Cappelli makes no attempt to link the parameters in his superpotential 
with the tricritical Ising Hamiltonian. 

A final issue with Cappelli's analysis is that he does not state along which curves one 
would obtain the (±) stable boundary states, he only says that they should be associated 
with the asymptotic regions (h b = ±00, 6 b = 00). The trouble is with his superpotential 
one would need to actually specify along which curves one should reach these asymptotic 
regions in. For instance if one reaches these limits along curves of the form h b = ±<5^ 
with < a < 1 then the stationary equation becomes 

+ 2S b m + m 3 = (4.6) 

dividing through by S b one gets 

3 

T^" 1 + 2m + ^ = (4.7) 

in the limit 5 b — > 00 and with < a < 1 the unique solution is m = which is the value 
associated to the (0) and not the (±) states. 



5 Conclusion 

In conclusion by extending the existing techniques of mean field theory we have presented 
a qualititive yet comprehensive picture of the boundary phase diagram of both Ising 
and tricritical Ising lattice models. Both showing agreement with the known pictures 
from conformal field theory. Furthermore we have found disagreement with the work of 
Cappelli et al. The crux of which is due to our equations are derived from first principles 
from the underylying lattice Hamiltonians whereas their work is based on what could 
apparentley be only an coincidental agreement between the nested pattern of RG flows 
in the boundary minimal models and the nested nature of the singular points in Arnold's 
theory of singularities. In taking this route, Cappelli et al. we think miss the true nature 
of the boundary critical points and go against the original principle of Landau theory, 
which is to only expand about small values of the order parameter. Finally the full mean 
field equations give useful information on the level probabilities throughout the phase 



27 



diagram, which we found useful in defining the stable (0) and (±) regions and in these 
regions mean field theory becomes exact. 

One could extend our work to all higher critical Ising models that would model the 
behaviour of their corresponding minimal models to give a qualititive agreement with 
the work of K.Graham [9]. Further as mentioned in the introduction that the principle 
motivation for the development of these techniques is to extend them to the case of 
conformal defects in minimal models as investigated in [10] and and to find nature 
of new defects. 
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